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Abstract 

Gravity modulation of an unbounded fluid layer with 
surface tension variations along its free surface is 
investigated. In parameter space of (wavenumber, 
Marangoni number) modulation has a destabilizing 
effect on the unmodulated neutral stability curve for 
large Prandtl number, Pr, and small modulation 
frequency, ft, while a stabilizing effect is observed for 
small Pr and large ft. As ft -> oo , the modulated neutral 
stability curves approach the unmodulated neutral 
stability curve. At certain values of Pr and ft multiple 
minima are observed and the neutral stability curves 
become highly distorted. Closed regions of 
subharmonic instability are also observed. Alternating 
regions of synchronous and subharmonic instability 
separated by very thin stable regions are observed in 
(l/ft^) space for the singly diffusive cases. 
Quasiperiodic behavior in addition to the synchronous 
and subharmonic responses are observed for the case of 
a double diffusive fluid layer. Minimum acceleration 
amplitudes were observed to closely correspond with a 
subharmonic response, X- m = ft / 2 . 


1. Introduction 

Under certain conditions such as thin liquid 
films or in a low-gravity environment, surface tension 
variations along a free surface may induce convection, 
referred to as theimocapillary convection. Orbiting 
spacecraft such as the NASA space shuttle and Russian 
Mir space station have proven to be useful 
environments for studying thermocapillary phenomena, 
which are often masked by gravity under terrestrial 
conditions. Thermocapillary phenomena are important 


to potential microgravity technologies such as crystal 
growing and materials processing applications as well 
as terrestrial applications such as coatings and drying 
processes. Onset of thermocapillary convection in the 
form of an extended Marangoni-Benard problem 
remains an active topic of study, in part due to its 
relative simplicity in terms of a motionless basic state 
and constant (or time periodic). The popularity of such 
models also stems from the understanding and insight 
they have provided to a broad range of physical 
phenomena (Pimputkar and Ostrach, 1981, Ostrach 
1982, Davis 1987, Legro et. al, 1990, and Kochmeider 
1993). Onboard space platforms the effects of vibration 
or g-jitter on thermocapillary systems is of particular 
concern (Nelson, 1991), although parametric excitation 
may also be of interest for terrestrial processes. With 
the question of g-jitter in mind, we focus on perhaps the 
simplest class of thermocapillary problems, the 
Marangoni-Benard problem, and consider how a fluid 
layer responds to time periodic accelerations or gravity 
modulation, which are imposed in the direction of the 
basic temperature gradient. 

Gershuni and Zhukhovitskii (1963) and 
Venezian (1969) considered the effect of temperature 
modulation on the Rayleigh-Benard problems where 
buoyancy drives convection. Two papers, Gershuni, 
Zhukhovitskii, and Iurkov (1970) and Gresho and Sani 
(1970), address gravity modulation of the Rayleigh- 
Benard problem. Gershuni, Zhukhovitskii, and Iurkov 
(1970) considered the linear stability of an unbounded 
fluid layer and fluid with in a vertical circular cylinder 
with modulation of both the mean vertical temperature 
gradient and the vertical component of gravity, while 
Gresho and Sani (1970) studied both linear and 
nonlinear behavior of the gravity modulated unbounded 
layer. In both studies, a one-term Galerkin expansion 
reduced the equations to a single second order ordinary 
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differential equation, ODE. Applying a direct analogy 
to the viscously damped pendulum with an oscillating 
vertical support, Gresho and Sani (1970) recast the 
problem in terms of the Mathieu equation reducing die 
number of parameters affecting the system behavior 
from six (fluid transport) parameters to three (Mathieu) 
parameters. Murray et. al. (1991) and Wheeler et. al. 
(1991) treated the gravity modulated onset of 
convection problem with respect to directional 
solidification. In these papers it was shown that certain 
combinations of modulation amplitudes and modulation 
frequencies, stabilized the unmodulated fluid layer that 
would normally be unstable, or destabilized a fluid 
layer whose unmodulated state is stable. 

More recently, Or and Kelly (1995) and Kelly 
and Or (1998) have investigated the effect of shear and 
temperature modulation on the Marangoni-Benard 
problem and noted interesting behavior that was also 
associated with surface deformation. For example, 
shear modulation by oscillating the lower (rigid) 
boundary destabilized the long wavelength mode but 
stabilized the finite wavelength modes. They found 
shear modulation could stabilize both finite and long 
wavelength modes while thermal modulation is 
effective in stabilizing the finite mode (Kelly and Or 
1998). 

The onset of convection (Rayleigh-Benard 
type) problem has been extended to double diffusive 
systems by Saunders et al. (1992) and Terrones and 
Chen (1993). In these studies density gradients due to 
temperature and species concentration gradients drive 
convection. Terrones and Chen (1993) studied the 
linear stability of gravity modulation and cross- 
diffusion on the onset of convection in an unbounded 
doubly diffusive fluid layer. They considered the Soret 
and Dufour effects as well as cases where the cross- 
diffusion effects were directly imposed. The results of 
both Saunders et. al. (1992) and Terrones and Chen 
(1993) demonstrate a rich variety of behavior, such as 
synchronous, subharmonic, quasi-periodic, and 
multiple minima. 

In this paper, parametric excitation of an 
unbounded fluid layer with surface tension variation 
along the free surface is considered and Floquet theory 
is applied to examine the stability of the modulated 
system. We initially focus on the gravity modulated 
Marangoni-Benard problem and then present a set of 
stability results for the modulated double diffusive 
Marangoni-Benard problem. The problem is first 
formulated for an unbounded double diffusive fluid 
layer where both density and surface tension variations 
due to temperature and/or concentration gradients 
induce convection. 


2. Formulation of Equations and Boundary 
Conditions 


An unbounded double diffusive fluid layer 
with dimension 0 < x] <> d is considered. The 
governing equations, continuity, momentum, energy, 
and species equations for incompressible flow are given 
below in Eqs. (l)-(4). 


SU 


Sx 


f = 0 


dt 1 9x ; 


-ldp‘ p . . 3 2 U* 

— TV—ftG )S.3 +V ... 

Po p„ aXjdX: 


3T* „. 5T* „ 9 2 T* 

__ 4- {J — = D.. ; — 

dt Sxj SxjSXj 


dc 


T - dc * 


0c* 


St J SXj SXjSXj 


( 1 ) 

( 2 ) 

(3) 

(4) 


where 8 y = Kronecker's delta and i=l ,2,3 and * 
denotes dimensional quantities. The dependent 
variables U* , T* , and c* are the velocity, 
temperature, and species, respectively. The time 
dependent body force term is periodic, of the form 

g*(t) = g* + g* cos^QY) . Density, p, is computed 
using the Boussinesq approximation, 
p = p D |l - P,(t* - T 0 ')— p 2 (c* - c 0 * )j , where p, is the 

thermal expansion coefficient and J } 2 is the solutal 
expansion coefficient. Reference values chosen for the 
buoyancy terms are the temperature and concentration 
basic state values of the lower surface, T 0 * and c * . The 
kinematic viscosity, v, thermal diffusivity, D n , and 
mass diffusivity, D 22 , are assumed constant. 

Impenetrable and no-slip velocity conditions 
are applied, and constant temperature and concentration 
are imposed at lower surface, X 3 = 0 . 

At X 3 = 0 (Bottom) 

U* = 0 , T*(0) = T 0 \ c*(0) = cj, (for i = 1 , 2 , 3 ) ( 6 a-e) 

The upper surface at x] = d is flat and nondeformable 
leading to equation (7a) while tangential stress balances 
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are given by(7b,c). Heat and mass flux 
given by equations (7d,e). 

conditions are 

At xj = d (Top) 


u; =0 

(7a) 

fsu; 3u;) 3T* 3c* 

H . . = Y.r-r-Yzr-r 

^3Xj 3Xj J 3Xj 3xj 

(7a,b) 

-pc, D „fC-Q\ -DafW 

(7d,e) 


where: j = 1,2; and p and c p are the dynamic viscosity 
and specific heat, respectively, with constant values. 
Surface tension, o, is approximated as a linearized 
function of the T* and c* , 
a*a 0 -yj(T* -T*)-y 2 (c* -c^.The surface 

tension variation with temperature, y x , and the surface 
tension variation with concentration, y 2 , are defined 
as y , = -(do / 9T*) c . p . and y 2 = -(9o / <3c - ) r p . , 

respectively (Adamson 1982). Q* and J* are heat and 
mass fluxes to the environment at the free surface. 

The velocity, temperature, and concentration 
basic state profiles are U* = 0 , 

r( X :) = T'(0)-AT*-^-, c’(x;) = c'(0)-Ac'^- Where 
d d 

the difference quantities of the form Ay* are defined as 
Ay* = y*(0)-y‘(d) . Following Joseph (1976), 

Saunders et. al. (1992), and Terrones and Chen (1993), 
the equations are linearized and then 
nondimensionalized. Reference values used to 
nondimensionalize the resulting disturbance equations 

ared, , AT*, Ac* for length, velocity, 

d D n 

time, temperature, and concentration, respectively. By 
assuming solutions of the form 

(u(x s ,t),e(x iy 0,^,0) = 

( w(x 3 , t), <K *3 > 0, X( x 3 » t)) expCi(a,x, + )) 
for the perturbation variables, velocity, u, temperature, 
0, and concentration, c, the following disturbance 
equations of x 3 -momentum, energy, and species are 
obtained. 

(d 2 -a 2 )w = 2 

-Pra 2 g(t)(Ra<|> + Rs/^ 22 ©) + Pr(D 2 -a 2 ) 2 w 

4 = (d 2 -a 2 j4> + w (9) 


0=2?j 2 (D 2 -a 2 )© + w (10) 

where g(t) = g 0 +g,cosf2t, g„ = /g c ,g, = gj /g c , 
and n = n'D„ /d 2 . 


Rigid, conductive, and permeable conditions, 
are imposed on the disturbance velocity, temperature 
and concentration at x 3 = 0. A flat-nondeforming free 
surface yields equation (12a) while the tangential stress 
balance is given by equation (12b). Equations (12c,d) 
are the disturbance flux conditions at the free surface. 


At x 3 = 0 , w = 0 Dw = 0 $ - 0 0 = 0 (1 la-d) 

At x 3 = 1 

w = 0 D 2 w = -a 2 (Ma<j> + Ms/ V n ©) (12a,b) 

D<t> + Nu<|> = 0 D0 + Sh0 = 0 (12c,d) 


where D denotes 9 / 9x 3 and ’ denotes 9 w / 9t . The 
resulting dimensionless parameters from the above 
equations are: the Prandtl number, Pr = v / D u , 
Diffusivity ratio, V n =D 72 ID n , thermal Rayleigh 
number, Ra = g c P t AT*d 3 / D ,, v , solutal Rayleigh 
number Rs = g c P 2 Ac*d 3 / D^v, thermal Marangoni 
number, Ma = y , AT*d / D u |i , solutal Marangoni 
number, Ms = y 2 Ac *d / D 22 p , surface Nusselt number, 


Nu = h,d / k , surface Sherwood number, 

Sh = h 2 d/D 22 . J 

The disturbance equations are reduced to a set 
of N ordinary differential equations, ODE’s, using a 
spectral (Chebyshev) collocation scheme. Floquet 
analysis is applied to examine the stability of the 
system of ODE’s (Meirovitch 1970, and Joseph 1976). 
The monodromy matrix is computed by integrating the 
N set of ODE’s N times. Floquet multipliers, p p which 
are the eigenvalues of the monodromy matrix are then 
computed, and the characteristic exponents, X p which 
determine the stability of the system are related to the 


floquet multipliers as X } = — Ln(pj) . The characteristic 


exponent with the largest real part, X, determines the 
stability of the system. If Re(X) is positive, disturbances 
grow, if Re(X) is negative, disturbances decay. The 
imaginary part of X is multivalued and characterizes the 


* hj and h 2 are surface heat and mass transfer 
coefficients defined in Pearson 1958, and Skarda et. al. 
1998 (among several refs.), k is the thermal 
conductivity of the fluid layer. 
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system response. The response is synchronous when 
Im (K) = nfJ , and subharmonic when 

\m(\) = + > where n is an integer value. When 

the two frequencies are incommensurate, |S| * n for 

lm(X) = 8“ , the response is quasiperiodic (Joseph 

1976, Saunders et. al. 1992, and Terrones and Chen 
1993). 

3. Results For Modulated Singly Diffusive Fluid 
Layer 


The effect of modulation frequency on neutral 
stability is first examined in (a,Ma) space, typically 
used to study stability behavior of the unmodulated 
Marangoni-Benard problem. The set of neutral stability 
curves in Figs, la-f corresponds to Pr, g„ and Ra values 
of 1, 5 and 1000, respectively. Below each curve, the 
system is stable and “small” disturbances decay, and 
above the curve, the system is unstable and the 
disturbances will grow in time. The unmodulated 
neutral stability curve is shown in each figure as a 
reference point to directly compare and contrast the 
effect of modulation on the neutral stability boundaries. 


Fig. la reveals the existence of two small 
local minima (denoted as humps) near the bottom of the 
neutral stability for Q=7. Ma c is associated with hump 

3 in Fig. la and shifts to hump 2 for £2=10 indicated in 
Fig. lb. This suggests that a double minima exists 
along the synchronous curve for some Q between 7 
and 10. In Fig. lc where Q =14.5, hump 2 forms a 
narrow finger that extends below the unmodulated Mac 
of 80, thus having a destabilizing effect on the 
unmodulated neutral stability curve. Another local 
minimum, hump 1, is found near a=l in Fig. lc. 
Examination of Figs, ld-f indicates hump 1 continues 
to grow as hump 2 recedes with increasing Q, 
suggesting the presence of another double minima. 
Hump 2 eventually disappears with increasing 
modulation frequency. Between Q of 17 and 1 8, a 
subharmonic closed neutral stability branch forms. The 
subharmonic region of instability grows with increasing 
modulation frequency until reaching an Q of 
approximately 35. The subharmonic branch then begins 
to shrink and shifts to higher wavenumbers and lower 
Marangoni numbers. The subharmonic branch 
continues to shrink and eventually disappears as Q 
increases to a value of 100 where only the synchronous 
branch is evident. For Q > 100 the subharmonic loop 



0 1 2 3 4 5 6 


a 



0 1 2 3 4 5 6 


a 



0 1 2 3 4 5 6 


a 


Modulated Neutral Stability Curve 

Unmodulated Neutral Stability Curve 

Figure 1 Neutral Stability Curves for Pr=l., g Q =0, gj=5, Ra=1000 
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shifts out of the range of Ma values explored while 
modulation is observed to have a stabilizing effect on 
the synchronous branch. The stabilization effect 
however diminishes with further increase of the 
modulation frequency. The shapes of the neutral 
stability curves also change with modulation. 

The behavior of critical Marangoni numbers, 
Ma^ and critical wavenumbers, a c , with respect to 
modulation frequency are shown in Figs. 2 through 3 
for Pr values of .01, 7.1, 10 and 100. In view of the 
complex topology of the neutral stability curves, the 
critical values for Pr=l are presented separately in Figs. 
4 and 5. This level of complexity was not observed in 
the neutral stability plots at the other Pr values (for the 
range of parameter values we investigated) therefore 
we discuss these Pr values first. 


increases, gravity modulation has a stabilizing effect on 
Ma c which reaches a maximum at some finite values of 
Q. The stabilizing effect then decreases with a further 
increase in Q, and approaches the unmodulated critical 
value of Ma c =79.607 as Q — ► a> . Values of Ma c , Q, 
and a c where maximum stabilization occurs are given 
in Table 1. 


Table 1 Values of Q, Ma c , and a c 
for maximum stabilization 


Pr 

g. 

Q 


<*c 

.01 

5 

3.4 

125.77 

2.68 

7.10 

5 

88.0 

123.67 

1.54 

10.00 

5 

116.0 

115.20 

1.60 

100.00' 

5 

1000.0' 

85.09' 

1.92’ 


* Maximum Stabilization Occurs Beyond Q=1000 



Figure 2 Critical Ma c versus modulation frequency, Q, 
for varying Pr values. Ra = -1000, go = 0, g, = 5. 



Figure 3 Critical wavenumber, a c , versus modulation 

frequency, Q, for varying Pr values. 

Ra= -1000^ = 0, g, = 5. 


Gravity modulation for sufficiently small Q, is 
destabilizing for all Pr values investigated. As Q 


Critical and extremum values of Ma and a for 
Pr=l are shown in Figs. 4 and 5. Critical Marangoni 
numbers and wavenumbers of the synchronous branch 
are shown as functions of modulation frequency in 
Figs. 4a,b, while extremum values of Ma c and a c for 

the subharmonic branch are shown in Figs. 5a,b. 



synchronous branch. ForPr=l 



synchronous branch. For Pr=l 
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n 


Figure 5a Ma c versus Q for the 
subharmonic branch. For Pr=l 



Q 


Figure 5b a c versus Cl for the 
subhaimonic branch. ForPr=l 

The location and critical values corresponding 
to the two double minima observed in Fig. 4a are given 
in Table 2. Physically, the double minima suggest 
possible occurrence of mode switching where the 
cellular (or roll) pattern alternates between two cell or 
roll sizes dictated by the two different wavenumbers. 
Such a process would be easily observed given the 
large difference in the two wavenumbers associated 
with each of the double minima in Table 2. 


Table 2 Location and critical values 
for double minima 


£2 

Ma,. 

a cl 

a c2 

9.00 

124.0 

1.7 

3.1 

25.24 

185.7 

1.2 

3.5 


Figs. 5a,b reveal that the subharmonic branch 
forms at an Q slightly less than 17.5 and disappears just 
beyond an Cl of 52. The subharmonic instability 
extends to Ma values less than -250, where stable 
temperature gradient occurs for the unmodulated 
problem. 


3.2 Modulated Marangoni-Benard & Rayleigh- 
Benard Instabilities in (l/Q^) Space 

For the unmodulated problem, neutral stability 
curves in (a,Ma) space for the Marangoni-Benard 
problem are qualitatively similar to those of the 
Rayleigh-Benard. Unfortunately this correspondence of 
stability behavior between (a,Ma) space and (a,Ra) 
space does not hold for the modulated problem. This is 
because changing Ra simultaneously effects both the 
steady instability agency (buoyancy) and the 
modulation amplitude which is also part of the 
buoyancy term. In contrast, varying Ma only effects a 
steady instability agency - that of the surface tension 
variation. The result is that stability behavior observed 
in the (a,Ra) plane is not likely to apply to similar 
conditions in the (a,Ma) plane. A more appropriate 
space to compare the modulated Rayleigh-Benard and 
Marangoni stability behavior is that of (l/Q^) space. 

Three neutral stability maps in (l/£2,g,) space 
corresponding to different values of Ma, and go are 
shown in Figs. 6a,b,c. The values of a, Ra, and Pr are 
2, 1000, and 1, respectively, for these figures. For go= 

1, the Rayleigh number of 1000 in Fig. 6a, is 49.2% 
higher than the unmodulated neutral stability Rayleigh 
number at a=2. In Fig. 6b, the mean gravity 
acceleration level, go, and Ma are zero: therefore both 
the steady instability agencies of buoyancy and surface 
tension are zero. And in Fig. 6c, the Marangoni number 
is 1 1 8.77, which is 49.2% greater than the unmodulated 
neutral value of Ma at a=2. 

Results in Fig. 6a correspond to the modulated 
Rayleigh-Benard problem with a rigid-conductive 
lower surface and a free-insulated upper surface while 
the steady gravity acceleration level, go, is set to one 
which corresponds to terrestrial conditions. Fig. 6a 
reveals that a finite acceleration amplitude is required 
to stabilize the layer. Further increase in g L eventually 
destabilizes the layer. In general, alternating bands of 
synchronous and subharmonic instability separated by 
thin stable regions are observed in Fig 6a which is also 
characteristic of temperature modulation results of 
Gershuni and Zhukhovitskii (1963) and Kelly and Or 
(1998). Gershuni and Zhukhovitskii (1963) referred to 
the lower most region of instability as the fundamental 
instability region and showed that it existed for Ra 
values exceeding the unmodulated neutral Ra value*. 


§ The situation equally applies to the Marangoni-Benard 
problem where Ma values replaces Ra values. 
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i/n 

(a) Pr=1, g 0 =1, Ma=0 



1/Q 

(b) Pr=1, g 0 =0, Ma=0 



1/0 

(c) Pr=1, g 0 =O, Ma=1 18.77 

\ •" ' j Subharmonic Instability 
ii 1 I Ml Synchronous Instability 


Figures 6a, b, c. Stability boundaries in (1/Q,g 3 ) space 
for a=2 and Ra=1000. 

In Fig 6b, the steady acceleration level is 
go = 0, corresponding to a zero gravity condition while 
all other parameter values and boundary conditions are 


identical those applied in Fig. 6a. Since mean instability 
agencies are zero for these results, i.e., Ma=go=0, no 
fundamental instability exists. The system is stable for 
all frequencies until a finite amplitude value is reached. 
The stable regions separating the fingers of instability 
become quite thin with decreasing frequency and 
increasing modulation amplitude. A similar set of 
results, but with free-conductive boundaries was 
presented by Saunders et al. (1992). They chose Ra 
and Pr values of 1000 and .1, respectively, and 
a = 7 t / 2 . While their Pr and wavenumber differ from 
the Pr and a values used for Fig. 6, the shape and 
orientation of the instability fingers are quite similar. 
Their minimum modulation amplitude for instability 
was 7.3 corresponding to an Cl of 10.5 as compared to 
the minimum amplitude of 4.05 and Cl of 24.4 for Fig. 
6b.” 

The gravity modulated Marangoni instability 
shown in Figs. 6c, exhibits a fundamental instability 
band for small values of g,. This is consistent with the 
Rayleigh-Benard instability discussed for Fig. 6a, since 
Ma exceeds its corresponding unmodulated neutral 
stability value, Ma=79.6. A minimum finite amplitude 
of 4.97 is reached at £1=20.0 where modulation 
stabilizes the fundamental instability. At larger 
amplitude, the resonant instability regions or fingers, 
are attained as in Fig. 6b. The appearance of the 
synchronous and subharmonic instability regions is 
qualitatively different from the modulated Rayleigh- 
Benard problem observed in Fig. 6a. In Fig. 6c, 
multiple regions of stability exist that engulf each of the 
subharmonic instability fingers, and a single region of 
synchronous instability, the fundamental instability, 
surrounds the multiple stable regions. In all plots in 
Fig. 6a,b,c, regions of alternating 
subharmonic/synchronous instability occur and are 
separated by regions of stability. 

4.0 Results For The Double Diffusive Fluid Layer 

The large parameter space associated with the 
gravity-modulated doubly-diffbsive Marangoni 
instability problem poses a serious challenge to a full 
characterization of such systems. For our doubly 
diffusive gravity-modulated problem there are twelve 
parameters (Pr, Ra, Rs, Ma, Ms, a, Nu, Sh, D^, go, g„ 


” For consistency, Saunders et al.’s (1992) 
dimensionless frequency was converted to £1 defined 
herein by Cl = Pr Cl Smndcn . 
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£2). n The difficulty in handling large parameter spaces 
was similarly acknowledged in the gravity modulated 
Rayleigh-Benard investigations of Saunders et. al. 
(1992) and Terrones and Chen (1993). For this study, 
we choose a single set of values for (Pr, Ra, Rs, 

Ms, Nu, Sh, gj, and again investigate stability behavior 
in the parameter spaces of (a,Ma) and (l/ft,g,). 

4.1 Stability Boundaries For Unmodulated Double 
Diffusive System 

For the unmodulated double diffusive 
problem, oscillatory (diffusive) neutral stability 
associated with a Hopf bifurcation is possible, in 
addition to the stationary (fingering) instability 
resulting from an exchange-of-stabilities that typically 
occurs for the singly unmodulated Benard problem. 
Modulated behavior also differs depending on the 
nature of the unmodulated instability, oscillatory or 
stationary. Therefore, we first compute stability 
boundaries for the unmodulated double diffusive 
problem using a normal mode analysis and use these 
results to choose parameter values that yield the desired 
behavior i.e. stationary or oscillatory. 

Stability boundaries are shown in (Ms^MaJ 
space in Fig. 7 for the parameter set (Pr,^ 22 ,Ra,Rs,Nu, 
Sh) with values (10, 0.1, 500, -500,0,0). The Pr value is 
representative of water-alcohol or thermohaline 
systems, although, for such systems is normally in 
the range, 10‘ 3 to 10' 2 . The shape of the stability 
boundaries is also characteristic of double diffusive 
systems with the small values (Legros 1990, 
Terrones & Chen 1993, Skarda et. al. 1998). A co- 
dimension two point, represented by the intersection of 
the oscillatory and stationary boundaries occurs at the 
point, (Ms^MaJ = (-9.0, 88.6). To the left of the co- 
dimension two point, the transition from stable to 
unstable behavior occurs via a Hopf bifurcation and 
oscillatory neutral stability behavior is observed. To the 
right of the co-dimension two point exchange of 
stability is observed resulting in the presence of 
stationary neutral stability boundary. Frequencies, X im , 
and critical wavenumbers, a c , associated with the 
oscillatory stability boundary are shown in the insert of 
Fig. 7. (Here, X is the eigenvalue with the largest real 
part obtained from a normal mode analysis of the 
unmodulated problem (Pearson 1958, Joseph 1976, 
Skarda et. al. 1998)). Both X im and a c increase with 


+t Even the singly-diffusive gravity-modulated problem 
required consideration of eight parameters. 


increasing Ma* values, moving to the left of and away 
from the co-dimension two point. 



Figure 7 Stability Envelope for Pr=10, 1^=0. 1 
(Insert shows variation of frequency, X m , and wavenumber, aj 

4.2 Gravity Modulated Double Diffusive Layer 

For the parameter values we chose that 
correspond to the stationary (fingering) instability, 

X m = 0 , the modulated stability behavior is similar to 
the singly diffusive results discussed in the previous 
section. However, for parameter values corresponding 
to oscillatory (diffusive) behavior, X m * 0 , of the 
unmodulated system, a rich variety of behavior is 
possible along the neutral stability curve. Only results 
for parameter values corresponding to oscillatory 
behavior are presented below. 

The effect of frequency modulation on neutral stability 
curves is shown in Figs. 8a, and 8b. The elements of the 
parameter set (Pr, *>22, Ms, Ra, Rs, Nu, Sh, go) are 
chosen as (10, 0.1, -700, 1000, -1000, 0, 0, 0), which 
correspond to oscillatory neutral stability behavior of 
the unmodulated problem. The unmodulated double 
diffusive neutral stability curve is provided in Figs. 8 
for comparison with modulated curves, and is 
represented by the dotted line. The corresponding 
frequency, X lim , values are shown in Figs.9a and 9b. The 
locations of the bifurcation points along the modulated 
neutral stability curves are tabulated in Table 3. 
Quasiperiodic behavior occurs when X im and Q are 
incommensurate. For the oscillatory region of the 
unmodulated layer, oscillatory behavior occurs, while a 
complicated temporal flow structure occurs for the 
quasiperiodic region of the modulated fluid layer in 
Figs. 8a and 8b. 
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Figure 8. Neutral Stability curves for gravity modulated doubly diffusion. 
Ms = -700, Pr = 10, V 22 = 0.1 , Ra = 1000, Rs = -1000, g 0 = 0, g, = 1. 

(a) 0 = 20, (b) 0 = 5. 




Figure 9 Imaginary part of most unstable characteristic root, X jm . 
Parameters, Ms, Pr, V 22 , Ra, Rs, g 0 , and g v have same values 
as those used for figure 8. (a) O = 20, (b) 0 = 5. 


Table 3 Location of Bifurcation Points For Neutral 
Stability Curves in (a,Ma) space Pr=10 £> 22 =.l 
Ra=500 Ms=-700 Rs=-500 go =0 g,=l 


n 

Bifurc Type 

a fc>ifr 

Ma** 

^im bifr 

20 

qu — > sub 

2.76 

122.5 

10.0 


sub qu 

4.48 

206.8 

10.0 


qu sub 

0.86 

208.0 

2.5 


sub -» qu 

1.61 

129.0 

2.5 

5 

qu-syn-qu 

2.18 

125.7 

5.0 


qu -» sub 

2.74 

126.1 

7.5 


sub -> qu 

2.98 

129.5 

7.5 


qu-syn-qu 

3.71 

157.0 

10.0 


* Subharmonic loop also observed for Q=50 


Two bifurcation points are observed along the 
doubly diffusive stability curve for Q-20 shown in Fig. 
8a. Moving from left to right along the neutral stability 
curve, quasiperiodic behavior gives way to 
subharmonic behavior which in turn bifurcates to 


quasiperiodic behavior again. Values of a and Ma 
corresponding to the bifurcation points are given in 
Table 3. Yhe critical Marangoni number, Ma^, occurs in 
the subharmonic region. However, careful examination 
of Fig. 8a also reveals the existence of a local 
quasiperiodic minimum. 

Reducing the modulation frequency further to 
Q=5 increases the complexity of behavior along the 
neutral stability curve. Six bifurcation points are shown 
in Fig. 8b. Five regions of quasiperiodic behavior are 
separated by intervals or points of subharmonic, 
synchronous, subharmonic, and synchronous behavior, 
respectively. The critical Marangoni number, Ma c is 
subharmonic, with local minima occurring in 
quasiperiodic and other subharmonic regions as well. 
Values of a, Ma, and X im , corresponding to the 
bifurcation points are given in Table 3. 

The frequency, A im , plots in Figs. 9a and 9b 
confirm that resonance is responsible for the 
bifurcations observed in Figs. 8a and 8b. Resonant 
behavior is expected when A im and Q are 


commensurate, i.e., X m 


n Q 


, where n is an integer 
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value. The maximum value of over the range 
0.5 < a < 5 for the unmodulated layer is 14.06, 
therefore, resonance is anticipated at some location 
along the neutral stability curve where the modulation 
frequency, Q, satisfies Q <28.12 . Subharmonic 
excitation is clearly visible in Figs. 9a and 9b as 


Q 

approaches the — value of 10 and 2.5 respectively. 


The smallest wavenumber, a, where X tm =5 in Fig. 9a is 
2.76, which is a quasiperiodic to subharmonic 
bifurcation point given in Table 3. Similar behavior 
occurs in Fig. 9b for ^=2.5. For Q= 5, both 
subharmonic and synchronous resonances occur for a 
values where the unmodulated X im of 2.5, 5, 7.5, 10 are 
obtained. 

The doubly diffusive Marangoni instability 
was also investigated in (1/Q,gi) space for a set of 
parameter values where unmodulated oscillatory 
behavior occurs. Values of a and Ma were chosen as 
2.2 and 111, respectively, while other parameters 
maintain their values prescribed above. This parameter 
set yields an oscillatory and slightly stable behavior in 
the absence of modulation where the eigenvalue with 
the largest growth rate was computed as 
X = - 0.2422 + i 6.137. 

The stability boundaries in (1/0, gi) space are 
shown in Fig. 10. A single continuous region of 
subharmonic instability was observed for sufficiently 
large modulation amplitudes, g P The subharmonic 
instability region extends to small g, values on either 
side of a small closed loop region of synchronous 
instability. A minimum amplitude, g, =0.228 is 
necessary to destabilize the layer. The parametric 

response, , in this case is subharmonic, where the 

corresponding modulation frequency is Q=12.3 which 
is twice the value of the unmodulated \ m indicated 
above. Quasiperiodic behavior was observed for 
Q>19.05. 


Time histories for the velocity, temperature, 
and species disturbances of a double diffusive fluid 
layer are shown in Figs. 11-13. Parameter values for all 
three cases correspond to neutrally stable systems. 
Oscillatory behavior for an unmodulated fluid layer is 
shown in Figs. 1 1 The subharmonic and quasiperiodic 
behavior shown in Figs. 12 and 13 correspond to 
parameter values chosen at two different locations 
along the neutral stability curve in Fig. 8b. 



Figure 10. Modulated double-diffusive 
Marangoni-Benard stability boundaries in 
(g^l/Q) space, a = 2.2, g 0 = 0, Pr = 10, 

= 0.1, Ra = 500, Rs = -500, 

Ma = 111, Ms = -700. 

5. Conclusions 

The stability behavior for both the singly and 
doubly diffusive Marangoni-Benard problems is more 
complex in the presence of modulation. For the singly 
diffusive problem, modulation had a stabilizing effect 
on neutral stability curves in (a,Ma) space at large 
modulation frequencies and approached the 
unmodulated neutral stability curve as Q -> oo . For a 
Pr of 1, both synchronous and subharmonic regions of 
instability were observed for certain values of 
modulation frequency. The synchronous branch 
became distorted with local minima and narrow fingers 
forming at certain values of Q. Sometimes both local 
minima were above the unmodulated Ma c , thus 

resulting in a stabilizing effect. At other times one local 
minimum extended below the unmodulated Ma^ 

Maximum stabilization, (largest MaJ occurred at some 
finite value of modulation frequency. 

In previous modulated Rayleigh-Benard 
studies, alternating regions of harmonic and 
subharmonic were observed in (1/Q,g,), however the 
shape of the unstable regions from those of the present 
modulated Marangoni study qualitatively differed. For 
the Marangoni-Benard problem, multiple regions of 
stability exist that engulf each of the subharmonic 
instability fingers, and a single region of synchronous 
instability, the fundamental instability, surrounds the 
multiple stable regions. 
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Figure 1 1 unmodulated oscillatory behavior Figure 12 subharmonic behavior Figure 13 quasiperiodic behavior 

Figures 11-13. Time histories of velocity, temperature, and species disturbances at the center of the layer, z=0.5, 
for neutral stability. Pr = 10, = 0.1, Ms = -700, Ra = 1000, Rs = -1000, & = 0, Q = 5. Figure 1 1, a = 2.5, 

Ma = 1 19.86. Figure 12, a = 1.2, Ma = 99.2777., Figure 13, a = 2.7, Ma = 125.8057. 


The gravity modulated doubly diffusive 
Marangoni instability was considered for a single set of 
Pr andZ^ values. For parameter values corresponding 
to the stationary (fingering) instability, the modulated 
behavior is similar to the singly diffusive results. When 
parameter values correspond to oscillatory (diffusive) 
behavior, X m * 0 , of the unmodulated system, regions 
of quasiperiodic, subharmonic and synchronous 
behavior were found to exist along the neutral stability 
curves. This behavior is consistent with that of the 
double diffusive Rayleigh-Benard instability reported 
by Saunders et. al. (1992) and Terrones and Chen 
(1993). As modulation frequency, Q, was decreased, 
bifurcations along the neutral stability curve became 
more frequent due to parametric resonances. For a 
given value of A. im , accompanied by a decrease in £2, 


clearly the number of bifurcations increases as n = 

(n is an integer value only). In (1/Q,gi) space, the 
minimum g, necessary to destabilize the layer was the 
result of subharmonic excitation and thus occurred at a 
modulation frequency twice the value of the 
unmodulated X ^ . 


References 

Adamson, A.W. 1982 Physical Chemistry of Surfaces, 
4 th edn. John Wiley. 


11 

American Institute of Aeronautics and Astronautics 











Davis, S.H. 1987, “Thermocapillary Instabilities,” 
Annual Review of Fluid Mechanics, Vol. 19, pp. 403- 
435. 

Gershuni, G.Z. and E.M. Zhukhovitskii, 1963, "On 
Parametric Excitation of Convective Instability,” 

PMM., Vol. 27, No. 5, pp. 779-783. 

Gershuni, G.Z., E.M. Zhukhovitskii, and I.S. Iurkov, 
1970, "On Convective Stability In the Presence of 
Periodically Varying Parameter," Journal of Applied 
Math. Mech., Vol. 34, pp. 442-452. 

Gresho, P.M., and R.L. Sani, 1970, "The Effects of 
Gravity Modulation on the Stability of a Heated Fluid 
Layer," Journal of Fluid Mechanics, Vol. 40, No. 4, pp. 
783-806. 

Joseph, D.D., 1976, Stability of Fluid Motions , 

Springer- Verlag, New York. 

Kelly, R.E., and Or, A.C., 1998, “ Influence of 
Temperature Modulation Upon The Onset of 
Thermocapillary Convection,” 36 th Aerospace Sciences 
Meeting, Reno, NV, AIAA 98-0652. 

Koschmieder, E.L., 1993, Benard Cells and Taylor 
Vorticies , Cambridge University Press. 

Legros, J.C., Dupont, O., Queeckers, P. Van 
Vaerenberg, S., 1990, “Double-Diffusive Convection 
And Its Effect Under Reduced Gravity,” In Low 
Gravity Fluid Dynamics and Transport Phenomena , ed. 
Koster, J.N. and R.L. Sani, American Institute of 
Aeronautics and Astronautics, Washington D.C. 

Meirovitch, L., 1970, Methods of Analytical Dynamics , 
McGraw-Hill Book Co. N.Y. 

Murry, B.T., S.R. Coriell, and G.B. Mcfadden, 1991, 
“The Effect of Gravity Modulation On Solutal 
Convection During Directional Solidification,” Journal 
of Crystal Growth, Vol. 1 10, pp.713-723. 

Nelson, E.S., 1991, "An Examination of Anticipation of 
g-jitter on Space Station and Its Effects on Material 
Processes, NASA TM 103775. 

Or, A.C. and Kelly, R.E., 1995, “Onset of Marangoni 
Convection in a Layer of Fluid Modulated by Weak 
Nonpianar Oscillatory Shear,” International Journal of 
Heat and Mass Transfer, Vol. 38, No. 12, pp^ 2269- 
2279. 


Ostrach, S., "Low-Gravity Fluid Flows," Annual 
Review of Fluid Mechanics, Vol. 14, 1982, pp. 313- 
345. 

Pearson, J. R. A., 1958, "On Convection Cells 
Induced by Surface Tension," Journal of Fluid 
Mechanics, Vol. 4, pp. 489-500. 

Pimputkar, S.M., and S. Ostrach, 1981, "Convection 
Effects In Crystal Grown From the Melt," Journal of 
Crystal Growth, Vol. 55, pp. 614-646. 

Saunders, B.V., B.T. Murray, G.B. McFadden, S.R. 
Coriell, and A. A. Wheeler, 1992, "The effect of Gravity 
Modulation on Thermosolutal Convection in an Infinite 
Layer of Fluid", Physics of Fluids A, Vol. 4, No. 6. 

Skarda , J.R.L., Jacqmin, D., M c Caughan, F.E., “Exact 
and Approximate Solutions to the Double-Diffusive 
Marangoni-Benard Problem With Cross-Diffusive 
Terms,” Journal of Fluid Mechanics, Accepted For 
Publication, 1998. 

Terrones, G., and C.F. Chen, 1993, "Convective 
Stability of Gravity-Modulated Doubly Cross-Diffusive 
Fluid Layers," Journal Fluid Mechanics, Vol. 255, 
pp. 301-321. 

Venezian, G., 1969, “Effect of Modulation On The 
Onset of Thermal Convection,” Journal Of Fluid 
Mechanics, Vol. 35, part 2, pp. 243-354. 


12 

American Institute of Aeronautics and Astronautics 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202*4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188). Washington, DC 20503. 


1. AGENCY USE ONLY ( Leave blank) 


TITLE AND SUBTITLE 


2. REPORT DATE 

June 1998 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Convective Instability of a Gravity Modulated Fluid Layer With 
Surface Tension Variation 


6. AUTHOR(S) 


J. Raymond Lee Skarda 


WU -962-24-00-00 


7. PERFORMING ORGANIZATION NAME<S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-11219 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM— 1998-207941 
AIAA-98-2599 


11. SUPPLEMENTARY NOTES 


Prepared for the 2nd Theoretical Fluid Mechanics Meeting sponsored by the American Institute of Aeronautics and 
Astronautics, Albuquerque, New Mexico, June 15-18, 1998. Responsible person, J. Raymond Lee Skarda, organization 
code 6712, (216) 433-8728. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Categories:34 and 29 


12b. DISTRIBUTION CODE 


Distribution: Nonstandard 


This publication is available from the NASA Center for AeroSpace Information, (301) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

Gravity modulation of an unbounded fluid layer with surface tension variations along its free surface is investigated. In 
parameter space of (wavenumber, Marangoni number) modulation has a destabilizing effect on the unmodulated neutral 
stability curve for large Prandtl number, Pr, and small modulation frequency, Q, while a stabilizing effect is observed for 
small Pr and large Q. As £2 — » oo, the modulated neutral stability curves approach the unmodulated neutral stability 
curve. At certain values of Pr and Q multiple minima are observed and the neutral stability curves become highly 
distorted. Closed regions of subharmonic instability are also observed. Alternating regions of synchronous and 
subharmonic instability separated by very thin stable regions are observed in (l/Q^) space for the singly diffusive 
cases. Quasiperiodic behavior in addition to the synchronous and subharmonic responses are observed for the case of a 
double diffusive fluid layer. Minimum acceleration amplitudes were observed to closely correspond with a subharmonic 
response, X im = D (2 . 


14. SUBJECT TERMS 

Microgravity; Marangoni; Benard; Onset of convection; Gravity modulation; g-jitter; 
Floquet; Stability; Thermocapillary; Hydrodynamic stability 

Bfl 

I _____ ___ __ 

17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 

18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 

19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 

20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 






















